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RESEARCH MEMORANDUM 


LONGITUDINAL FREQUENCY-^SPONSE CHARACTERISTICS OF A 
35° SWEPT-WING AIRPLANE AS DETERMINED FRCM FLIGHT 
MEASUREMENTS, INCLUDING A METHOD FOR THE 

EVALUATION OF TRANSFER FUNCTIONS 
By William C. Triplett and G. Allan Smith 

SUMMARY 

Longitudinal frequency— response characteristics of a 35° swept-wing 
airplane computed from flight measurements are presented herein. 
Pitching— velocity and normal acceleration responses to pulse— type eleva- 
tor inputs have been analyzed to obtain frequency responses and numeri- 
cal values of the coefficients of the transfer functions. These coeffi- 
cients are plotted as functions of Mach number over the Mach number 
range of 0.59 to 1.05. All flight data were recorded at an altitude of 
approximately 35^000 feet, and effects of altitude variations on the 
frequency— response characteristics are calculated. 

Also presented are the pitching— velocity frequency responses to 
displacements of the adjustable stabilizer. These were computed using 
experimentally determined values of the ratio of stabilizer effective- 
ness to elevator effectiveness. 

A graphical method for determining transfer coefficients from 
frequency— response data is described. This method eliminates lengthy 
numerical culculations usually associated with this type of problem. An 
electronic analogue computer was used in conjunction with the above 
method so that the graphically determined results could be modified 
until the output of the computer matched the original flight data. 


INTRODUCTION 


In the design of automatic— control equipmsnt for high performance 
aircraft, it is necessary to consider the dynamic— response 
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characteristics of the aircraft. These characteristics can he conven- 
iently expressed as transfer functions which are expressions that des- 
cribe the motion of the airplane for the various flight conditions of 
interest. Often these expressions can be predicted by using stability 
derivatives obtained from theoretical studies and from wind-tunnel tests. 
In many cases, however, particularly in the transonic speed range, it is 
necessary to resort to actual flight— test procedures in order to deter- 
mine the airplane frequency response from which the desired transfer 
functions may be obtained. 

Presented herein is one phase of a flight— test program being con- 
ducted by the NACA for the puirpose of determining longitudinal, lateral, 
and directional frequency-response characteristics of a 35° swept— wing 
airplane. This report considers only the longitudinal case and presents 
variations of the coefficients of the transfer functions over the Mach 
number range of 0,59 to 1.05 for one flight altitude. The results of 
this phase of the test program were determined from flight data similar 
to that described in reference 1 and this report may be considered as an 
extension thereof, 

A second objective is the presentation of a general method for 
obtaining transfer functions from frequency— response data utilizing a 
special set of templates and an electronic analogue computer. The tem- 
plates afford a neans of rapid determination of the form of the transfer 
functions as well as good approximations of the numerical coefficients 
involved. The analogue computer may -then be used to modify these mrnBr— 
ical values and to check them against the flight measurements. 


NOTATION 



lift coefficient 


C. 


m 


pitching-moment coefficient about airplane center of gravity 


I 


y 


moment of inertia about the lateral axis through the center of 
gravity, slug— feet squared 


El, Kz amplitude scale factors 


M 


Mach number 


S 


wing area, square feet 


T 


time interval, seconds 


Y 


flight velocity, feet per second 


a 


sonic velocity, feet per second 
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c 

g 

m 

n 

P 

t 

a 

6 

P 

2 

T 

(Ji 

“n 

a 


mean aerodynamic chord of wing, feet 
acceleration due to gravity, feet per second squared 
angle of incidence of horizontal tail, radians 
mass, slugs 

normal acceleration, feet per second squared (unless other- 
wise noted) 


the operator 


dt 


a 




pitching Telocity, radians per second 

time, seconds 

angle of attack, radians 

eleTator deflection, radians (unless otherwise noted) 
atmospheric density, slugs per cubic foot 

dynamic pressure, pounds per square foot 

damping ratio 

characteristic time, seconds 
frequency, radians per second 

undamped natural frequency, radians per second 

da 

dt 

dt 

dCL 

da 

d6 
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C™ 

da 

Cmg 

dCm 

dF" 

”q 

dCm 

d(qc/2Y) 


dCm 

d(dc/2V) 

Ma 


Md 

1 <=% 

Mq 

i pvsc^ 

Za 

- \ Ci^ 

Zb 

1 °Lb 


TEST EQUIPMENT 


The test airplane was a standard North Arerican F86A— 5 with exter- 
nal booms added as shown In figures 1 and 2. The physical character- 
istics and instrumentation of this airplane are described in reference 1. 
Of particular Importance to the analysis used in the present report is 
the dynamic behavior of the instruments used to record pitching velocity 
and no rmal acceleration. Pitching velocity was measured by a rate gyro 
with direct optical recording. The natural freQ.uency of this instrument 
was 10.0 cycles per second and the dancing ratio was O.67. Nomal accel- 
eration was measured by a standard NACA air— damped vane— type accelero- 
meter with a natural frequency of l4.7 cycles per second and a damping 
ratio of I.95. 
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TEST PEOCEDUEE AND ANALYSIS 


The procedure used In the determination of the airplane frequency 
response and corresponding transfer functions can he divided into four 
distinct steps; namely, flight measurements, frequency— response computa- 
tions, determination of transfer coefficients, and a final check utiliz- 
ing an analogue computer. These steps are outlined in the following 
sections. 


Flight Measurements 


Transient responses of pitching velocity and normal acceleration to 
pulse— type elevator inputs were recorded over the Mach numher range of 
0,59 to 1 . 05 , All flights were made at an altitude of 35^000 feet at 
values of trim lift coefficient that varied from O .35 to 0,12, Sanq)le 
time histories are shown in figure 3« 

To insxire consistent data, the airplane was trimmed so that varia- 
tions in speed and altitude during each rim were held to a Ttilninnun- To 
attain test speeds greater than 0,95 Mach numher, it was necessary to 
dive the airplane, thus causing small altitude changes during these runs. 
In all tests, after applying the elevator pulse, the control was held 
fixed until the transient response was completely damped. As a check on 
accuracy, at least two flight records were tsiken at each speed. 

The magnitudes of the applied elevator poilses generally varied from 
10° to 12°, measured from trim position as shown in figure 4. It is 
desirable to use large inputs so that the transient responses can he 
more accurately meas-ured. On the other hand, the analysis methods pre- 
sented herein are valid only for linear systems. The input magnitudes, 
therefor, should he limited to a range in which the pertinent aerody- 
namic derivatives are essentially linear. Figure 4 shows the trim ele- 
vator angle at each flight— test condition in addition to the absolute 
magnitude of the input piilse. Figure 5 shows the corresponding trim 
angles of attack and the maximum and minimum angles reached during the 
oscillations that followed each elevatcxr pulse. 

A transient— type input, in particular a pulse, was used in these 
tests mainly because it enables an entire airplane frequency response to 
he computed from one flight run of a few seconds duration. Sinusoidal 
elevator inputs could have been used to determine the steady— state fre- 
quency response (amplitudes and phase angles) directly and thus elimi- 
nate the calculations described in the following paragraphs. The use of 
this type of input was not considered to be feasible because of the 
large amount of flight time required inasmuch as only one point on the 
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frequency-response curve could le obtained from each flight rim. In 
Q'^^i'^ion it would require the installation of special sine^^^ave generat — 
ing equipment in the airplane. 


Frequency-Eesponse Computations 


The next step in the analysis is the conversion of the transient 
time histories into the frequency domain. 

Under certain conditions a time function such as the pitching 
velocity q(t) can be transformed into a complex frequency function 
l(ito) by means of the Fourier integral relation 

q(iO)) = r q(t) e~^“^ dt (l) 

Jq 

This integral must be eYaluated from time zero to infinity for each fre- 
quency (a 3 at which q(iw) is desired. Obviously, the integration can 
be accomplished only if the behavior of q(t) is known for an infinite 
time. Since q(t) can be measured only for a finite time T, it is 
necessary for the system to reach steady-state conditions before time T 
such that q(t) may be expressed analytically between the time limits T 
and infinity. In addition, the product q(t)e”^^^ must converge as t 
approaches infinity. 


Examination of the data used in this analysis showed that in every 
case the forcing function &(t) and the responses q(t) and n(t) 
became constant after a finite time interval T. These steady— state 
values may be expressed as &p, qip, and np, respectively, and thus 
equation (l) may be divided into two parts as follows: 



The second integral can be evaluated analytically so that 


5(lu) = S 



( 2 ) 


For computational purposes equation (2) may be divided into its real and 
Imaginary parts represented by the symbols E.P. and I.P., respectively, 
so that 


q(iu) = (R.P.) + i(l.P.) 
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where 


E.P. 


I.P. 


The two integrals involving the transient part of q(t) may he evalu- 
ated by any of several approxiinatlon methods. In analyzing the data of 
this report, Simpson's rule was used to find the area under the product 
ciirves q(t) cos ojt and q(t) sin o)t. Values of q(t) were tabulated 
at 0.05— second intervals and the integrations were carried out at each 
integral value of frequency from 1 to 10 radians per second. 

After obtaining the numerical values of E.P. and I.P. at each fre- 
quency, q(iw) can be expressed in polar notation such that 

q(iw) = |q| e^^q 

where the amplitude, |q| = ^Ce.P. + (I.P. and the angle, 

= arc tan I.P./R.P. 

The above integration process was repeated for the elevator forcing 
function 6(t) to determine 

6(ito) = |6|e^^^ 

Then the ratio of c[/5 was determined as 

I - llSI = I § I - 1 § I 

where |q/6 | is the amplitude ratio and 9- /g is the phase difference 
between the two quantities, 9^-96* ^ 

Similar calculations were made using the no rmal acceleration tine 
histories in order to obtain the qixantities |n/6| and 9n/6* 

The final result of the above calculations is a graphical repre- 
sentation of the airplane transfer function (frequency response) for a 
given test condition. This is plotted as curves of amplitude ratio and 
phase angle versus frequency. The remainder of this section of the 
report is concerned with the determination of analytical expressions for 
these curves. 


qm r / X 

= — sin wT+ / q(t; cos ut dt 


= — ^ cos coT — / q(t) sin mt dt 
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Eeference 2 is one of many sources that explain in more detail the 
use of the Fourier transform for prohlems of this type. 

Determination of Transfer Coefficients 


A graphical method was used to find the type of transfer functions 
that best define the measured q/6 and n/6 frequency responses. 
Numerical values of the transfer coefficients were also determined in 
the same operation. This method involved the use of a set of templates 
developed by C. S, Draper of the Instrumentation Laboratory of the Massa- 
chusetts Institute of Technology. These templates are described fully 
in reference 3. The design of the templates was based on the following 
considerations : 

1. A rational function of complex frequency p can, in general, 
be expressed as 

F(p) = K » » » (p+an) / x 

(p+bi)(p+bs) , . . (p+bj^) 


where the a^^ and b^ may be either real or complex, and, if complex, 
^■PP®^^ conjugate pairs and where K is any real number. 

Vfhen aj^ is a real number, a factor (p+a^) may be expressed as 

I (1+Tp) 

where t = ~. When a^ and a^^^^ are a complex conjugate pair the 

two factors (p+8-i)(p+aj^^j^) can be expanded, by introducing new coeffi- 
cients, into the form (p^+ 2 Cu^p+u)j ^2 j ^ a±eii+j_ and 

2Cun = aj + a^+i. Factoring out gives 


f 1 + 2C -^ + 

V tJn Wn^y 

The terms ^ and that appear outside the parentheses, being real 

numbers, may be grouped together with the multiplying factor K, When 
an aj or bi is zero then a factor p will be present in either the 

nu^rator or the denominator. Then, F(p) in factored form will consist 
Gnuiroly of* comb 1110,13 ions of* bh© bhi*©© fopins* 


Since we are interested only in the steady-state frequency response 
the coD^lex variable p can be replaced by the frequency variable iu, 
and each of the three factors may be written in polar form as follows: 
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For (l+Tp ) , 


1 + xio) = a/i+t^u^ 

where 


ih) 


cpi = tan~^ TO) 





and for p 


iw = we^'Ps 


(5) 


(6) 


where 


cpg = tan"^“ = 90° 

2. When to) is taken as a nondlmensional fregiiency variable all 
first-order terms of the form (l + Tico) can he defined hy a single 
amplitude curve and a single phase-angle curve as shown in figure 6. 

In a similar manner all possible second— order terms of the form 
2 

can be defined by a family of curves with co/cjj^ as 

the nondimensional frequency variable. There will be a different pair 
of curves for each value of damping ratio as shown in figure 7. The 
factor iw has not been plotted because (from equation (6)) its ampli- 
tude is simply equal to o) and the phase angle is a constant 90 ^. 


1 + 


Icj 


ICO 


3 . Since a typical transfer function may consist of mcjre than one 
of these factors in both numerator and denominator, it is of great advan- 
tage to plot the amplitudes to a logarithmic scale so that multiplication 
or division of factors of the type shown in equations (4) to (6) may be 
accomplished by mere graphical addition of the amplitudes. The nondi— 
mensional frequency has also been plotted to a logarithmic scale as 
shown in figures 6 and J, Then as frequency approaches either zero or 
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infinity the amplitudes become asymptotic to straight lines. In both 
figures as frequency decreases the amplitudes approach asymptotically 
the value of unity. At high frequencies the amplitude of (l + Tiu) is 
essentially equal to tw and thus it plots as a straight line of slope 1 
(fig. 6). Similarly, the amplitude of a second-order term approaches 
the value (u/Uq)^ as frequency increases and is asymptotic to a 
straight line of slope 2. Since a second— order term of this type nor- 
mally appears in the denominator, it has been plotted as 


and thus the slope of the amplitude curve is -2 (fig. 7). In both plots 
the asymptotes Intersect at unity on the nondimensional frequency scale. 
This Intersection is termed "the breakpoint." 

The phase angles are plotted to a linear scale because when factors 
of the types shown in equations (k) to (6) are combined, the resultant 
phase angle of F(p) la merely the algebraic sum of the individual 
angles. It can be seen in figures 6 and 7 that the angles of these 
first— and second-order terms approach 90° and — 180°, respectively, at 
high frequencies. Thus an an^ilitude slope of 1 corresponds to an angle 
of 90°, while a slope of -2 corresponds to an angle of — l80°. In gen- 
eral, as u approaches infinity the amplitude of F(p) approaches 

Bj(n-m), vhere n and m are defined in equation (3). Then the ampli- 
tude slope is equal to (n-m) on the logarithmic plot and the phase nnpl a 
I s 90(n-m)°. 

The templates used in the analysis of the data were accurate repre- 
sentations of the curves shown in figures 6 and 7 cut from transparent 
®9'^©rial with the breakpoints marked. Twelve pairs of second— order tem- 
plates were included in the set to give values of damping ratio from 
0.1 to 1.0 in increments of 0.1 and in addition values of 0.05 and 0.15. 

The computed frequency-response data (amplitude ratios and phase 
angles) were plotted to the same logarithmic scale as the templates, and 
then by a trial-and-error approach the template or combination of tem- 
plates that best matched the given data was determined. To determine 
the proper combination, various aiDplitude templates were positioned in 
such a way that their algebraic sum matched the experimental amplitude 
ratio cxirve. For each amplitude template the corresponding phase— angle 
template was placed on the experimental phase-angle curve so that the 
breakpoint of each pair was alined with respect to frequency. The 
various pairs of templates were then adjusted until the sum of the indi- 
vidual amplitudes and phase angles most nearly matched the test data. 
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noted and then the appropriate values of t and were determined 

directly, since for each first order term 


and for each second— order term 

n "Bp 

In general, each factor will have a different hreakpoint. 

At first glance, this procedure may seem extremely tedious hut it 
has been found that with practice one can determine coefficients quite 
rapidly in this manner. For a more complete description of the princi- 
ples involved in this type of graphical representation reference 4 is 
recommended. 

Final Check 



As a final step in the caj-culations a Reeves analogue computer was 
used to check the results of all previous computations. This was accom- 
plished hy placing the actual time histories of elevator-control motion 
into the computer hy means of an input table. Then this input was fed 
into a circuit representing the equations of motion as obtained in the 
previous step and the outputs of the machine were obtained. These out- 
puts were compared to the actual time histories of pitching velocity and 
normal acceleration obtained in flight. By changing dial settings on 
the computer the transfer coefficients could be adjusted until the out- 
put of the computer most nearly matched the actual flight data. (Such 
a comparison is shown in fig. 3») Thus, in addition to checking the cal- 
culations, this step actually resulted in a refinement of the coeffi- 
cients obtained with the templates. In general, operations of this type 
can be conveniently accomplished on the computer only when the form of 
the transfer function is known. 


RESULTS AND DISCUSSION 


The results presented in the following paragraphs were obtained 
from transient time histories such as shown in figure 3 ty the methods 
described in the previous section. Discussed here are the airplane 
frequency— response curves and the actual transfer functions that define 
these curves. Also presented as supplementary results are the Mach 
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nuaiber variations of elevator effectiveness and the effect of altitude 
on the transfer coefficients. 


Frequency— Response Curves 


Plotted in figures 8 and 9 are frequency— response amplitude and 
phase— angle curves of both pitching— velocity and normal— acceleration 
responses to elevator inputs for several different Mach numbers. 

Although two records were analyzed at each flight speed, the results of 
only one have been plotted. The agreement between the two results at 
each Mach num'ter was generally good with variations in amplitude ratio 
rarely exceeding 5 percent of their maximum values, while the differences 
in phase angle were usually less than 5°. 

The pitching-velocity frequency-reponse amplitudes of figure 8(a) 
show an increase in natural frequency with Mach number and also a reduc- 
tion in amplitude ratios above a Mach munber of O.76, indicating a 
decrease in elevator effectiveness. The normal— acceleration frequency- 
response amplitudes of figure 9(a) show, in general, the same character- 
istics with some irregularities at the higher flight speeds. 

Figures 8(b) and 9(t) show the corresponding phase-angle variations with 
Mach number. 

Each transient reponse was analyzed at integral values of frequency 
from 1 to 10 radians per second. In some cases, however, there appeared 
to be Inconsistencies in the calculated responses at either end of the 
frequency range where the magnitude of the response was small. This was 
particularly evident over the low-frequency parts of the normal- 
acceleration responses obtained at the higher flight speeds. For this 
reason, as well as for clarity, portions of some of the responses have 
been omitted from figures 8 and 9» In general, the pitching— velocity 
data are considered more reliable than the acceleration data because the 
dynamic characteristics possessed by the rate gyro were superior to those 
of the accelerometer. 

In figure 10 is plotted a series of amplitude— ratio curves of 
pitching-velocity responses to inputs of an adjustable stabilizer con- 
trol for various Mach numbers. These plots were obtained by multiplying 
the amplitude curves of figure 8(a) by values of d5/di+ obtained from 
static flight tests and shown in figure 11. Figure 10 illustrates the 
type of responses that would have been obtained if the stabilizer had 
been used as the maneuvering control and there is no indication of the 
large decrease in coxitrol effectiveness previously noted in figure 8(a). 
For these calculations the phase angles for the q_/U response are the 
same as those shown in figure 8 (b). 
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Transfer Functions 


The airplane transfer f\mctions were determined from q./5 and n/5 
frequency responses (figs. 8 and 9) l^y using the templates described in 
the Analysis section. Final values of the coefficients were obtained 
from the analogue ccaaputer by matching the original time histories. The 
transfer functions for this particular airplane were found to be of the 
following form; 




= El 

' 1 + 2C 

= Ks 

1 + 2C 


1 + Tiu 



1 




As indicated by figure 3^' these equations adequately define the tran- 
sient motion of the airplane as measured in flight. These equations sire 
also of the same form as the theoretical transfer functions developed in 
appendix A. 

The Mach number variations of the coefficients Ki and 

Ks are plotted in figures 12 and 13* The damping ratio curve of 
figure 12 was faired by referring to the damping coefficient curve 
plotted in reference 1 where more data were available to define the 
sharp variations between Mach numbers of 0.88 to 0.95. An indication of 
the consistency of the measurements and calculations is given by the 
scatter of the points at each flight speed. The coefficients ^ and 
(Jq were determined directly from the q/6 response and then were veri- 
fied using the n/6 response. 

In determining the transfer coefficients of the q/S response 
there was some uncertainty as to the best values for t and Ki. It 
was found that in some cases the test data could be equally well satis- 
fied by a range of values of these two coefficients. For examplej a 
decrease in t could be compensated for by an increase in Ki such that 
the product Kit remained constant. Examination of the q/£ transfer 
function will indicate the reason for this indeterminacy. When 
expressed in polar form the numerator of the amplitude ratio |o./8| is 

KiVl + t2o)2 


It can b© seen that as frequency increases becomes large compared 

to unity^ and the numerator may be approximated simply by Kitco. Thus 
it appears that KiT is the most essential quantity in the numerator of 
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the q./& transfer function and this has been plotted in figure l4. The 
resiilts in this form exhibit less scatter than the individual coeffi- 
cients Ki and t shown in figure I 3 . 


Elevator Effectiveness 


Also plotted in figure l4 is the variation with Mach number of the 
elevator effectiveness Cmg computed directly from the product Kx''" as 
described in appendix B. This quantity is compared to wind-tunnel 
values obtained from reference 5, and also to unpublished static fli^t 
data. The agreement between the three curves is generally good. There 
is scmie scatter, however, in the test points computed herein which may 
be attributed to the nonlinearities that are discussed in the section 
on accuracy. The method used here is not suggested as a neans for deter- 
mining Qnig since this can usually be obtained more directly from 
static flight measiirements. However, the computation is Justified as an 
additional check on the validity of the nethods used in this report. 


Altitude Variations 


Since the data presented in this report were all obtained at an 
altitude of approximately 35,000 feet, calculations, as descfibed in 
appendix C, were made to show the effects of changes in altitude on the 
freq^aency— response characteristics of the airplane. The results of 
these calculations are plotted in figure 15 , which shows the variation 
at constant Mach number and at constant lift coefficient of each of the 
transfer coefficients over an altitude range of from sea level to 
50,000 feet. Calculations based on these curves would be valid only at 
Mach numbers where altitude variations do not result in departures from 
the linear range of the lift curve. 

In figure I 6 are plotted typical pitching— velocity frequency- 
response curves calculated for the test airplane at a Mach number of 
0.81 for altitudes of sea level, 15,000 and 35,000 feet. It can be seen 
that as altitude increases both the amplitude and natural frequency are 
reduced. Calc\ilations for normal-acceleration response would show simi- 
lar tendencies with possibly even greater changes in amplitude. 


ACCURACY 


In order to assess the accuracy of experimentally determined trans- 
fer coefficients, consideration must be given to certain basic factors 
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that affect the solution of this type of problem. The templates vlll 
indicate the type of transfer function that best fits a set of experi- 
mental data over a particular frequency range. Changing the width of 
the frequency spectrum analyzed will, in many cases, not only change the 
values of the coefficients but also change the type of the transfer 
function. This would be true in general of a system that has more than 
one well-defined mode of motion. The airplane falls into this category. 

In the analysis of this report only a narrow frequency band was 
considered. This band bracketed the resonant peak of the short— period 
mode and thus, as would be expected by theory, the templates showed the 
system to have essentially two degrees of freedom (a second-order trans — 
fer function). If flight data had been obtained that would be suitable 
for analysis at lower frequencies, then we would have expected the 
system to be more closely represented by a fourth-order transfer func- 
tion that would have included the long— period (phugoid) mode in addition 
to the short— period mode mentioned above. Similarly, if the analysis 
had included very high frequencies, one or more vibration or flutter 
modes might have further complicated the resultant transfer function. 

Other factors that affect the acciaracy and consistency of the com- 
puted transfer coefficients are discussed in the following paragraphs. 


Nonlinearities 


The analysis used in this report is based on the assumption of 
linear behavior of all aerodynamic derivatives, thus implying that the 
frequency response and the related transfer coefficients are completely 
independent of the shape and magnitude of the forcing function. In cal- 
culating the dynamic^'esponse characteristics of aircraft, this assump- 
tion is generally valid, provided that only small disturbances from 
trimmed flight are considered. There may be cases, however, where non- 
linear derivatives exist even under these conditions. Then, at each 
fli^t speed and altitude, the computed frequency responses will vary in 
some manner with the type of control input used and also with the magni- 
tude of the input. This will lead to corresponding variations in the 
calculated transfer coefficients, and in extreme cases it may be impos- 
sible to find a set of linear- coefficients that adequately describe the 
transient flight data. 

Investigation has shown that, in the test airplane, the variation 
of pitching moment with elevator angle tends to become nonlinear 
(Cjug decreases) as the angle is increased above a certain value. 
Unpublished data have demonstrated this fact at a Mach number of 0.80 
where airplane responses to pulses of varying magnitude were recorded. 

It was found also that the amplitude ratio of the frequency response 
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decreased slightly as the magnitude of the pulse was Increased. The 
only coefficients appreciably affected were Ki and Ka. The varia- 
tions In Ki resulted In similar variations In the calculated values of 
Cmg which are in effect linear approximations of the true nonlinear 
curves. 

The results of this report indicate the same nonlinear effects. At 
several different speeds (such as O.69 Mach number) variations in the 
coefficients Ki and shown in figure 13» were directly traceable to 

variations in the sizes of the input pulses (fig. 4). For the most part 
these effects were small. 

From the equations developed in appendix A it can be seen that, 
through the frequency range of interest, both |q/5| and |n/6| are 
approximately proportional to 0^^. Thus it is reasonable to expect 
only the multiplying factors Ki °and to be affected by the simple 

nonlinearity in Cmg. 

There was no evidence of nonlinearities in any of the other stabil- 
ity derivatives. 


Instrumentation 


The dynamic characteristics of the recording InstruBBnts are of 
great importance. Each instrument has its own frequency response which 
must be considered in determining the true airplane response. Further- 
more, the natural frequency of the instrument definitely limits the 
highest frequency at which the data may be analyzed. 

Any tendency on the part of a recording instrument to drift with 
time would result in erroneous steady-state values of the airplane res- 
ponse (errors in the low-frequency region). 

The rate gyro used in this test program was found to have satis- 
factory djmamic amplitude characteristics, and oorrections were applied 
account for instrument lag. The linear accelerometer had less 
desirable dynamic characteristics because of its extremely high damping 
ratio (1.95) which normally means that amplitude as well as phase-angle 
corrections must be made^ but, because the normal acceleration data are 
of only secondary importance to this report, the performance of the 
instrument was considered to be adequate even thou^ amplitude correc- 
tions were not made. There was no evidence of drift with time in either 
InstruBBnt . 
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Flight Technique 


As mentioned previously, it was necessary to dive the airplane in 
order to obtain records at the highest speeds, which resulted in changes 
in both speed and altitude. These changes increase the difficulty in 
finding a set of coefficients that will describe the response of the 
airplane. Inadvertent disturbances such as those caused by rough air 
will cause airplane motions that are not consistent with the control 
input and thus may result in serious discrepancies in the airplane fre- 
quency response. 


Data Reduction 


The most likely source of error in the data reduction is the film- 
reading operation. However, such random errors should only affect the 
high-frequency part of the response, provided that the steady— state por- 
tion of the transient records is carefully measured. 

The use of Simpson *s Rule in obtaining the Fourier integrals of the 
transient records was found to be a highly accurate method, provided 
that the time history is evaluated at small enough time intervals. 
Unpublished work has shown that if a time interval of 0.05 second is 
used, then transient records of the type shown in this report can be 
analyzed to frequencies as high as 12 radians per second with errors of 
less than 1 percent. In general, the time interval should be small 
enough so that each cycle of the product curves q(t) sin wt and 
q(t) cos (jt can be defined by at least six measurements* 

In using the templates to determine transfer coefficients it was 
generally not possible to exactly match the damping ratio of an experi- 
mental frequency-response curve with any of the available second-order 
templates. This type of error could be minimized by expanding the set 
to include more intermediate values of In any event the damping 
ratio can be corrected by making use of the analogue computer as n©n- 
tioned previously. 


CONCLUDING REMARKS 


The primary objective of this report is the presentation of dynamic 
longitudinal— stability characteristics of a 35^ swept— wing airplane as 
determined from transient flight data. A secondary objective is the 
presentation of a method that enables these characteristics to be 
expressed in the form of airplane transfer functions for each fli^t 


CONFIDENTIAL 


18 


CONITDEKTIAL 


MCA EM A51G2T 


^ condition. Time histories of pitching— velocity and normal-acceleration 
responses to pulse-type elevator inputs were satisfactorily analyzed by 
means of the Fourier integral transformation to obtain the airplane fre- 
quency responses in graphical form. The equations or transfer functions 
that most nearly matched these cxirves were then readily determined by 
means of special templates. The analogue computor was useful in obtain- 
ing more satisfactory numerical coefficients of the transfer functions 
and in checking all intermediate calculations against the original 
flight measurements. 

E x a m i n ation of the pitching- velocity frequency-response curves for 
elevator-control inputs shows decreasing anq)litude ratios as Mach number 
is increased above 0,76, indicating low elevator effectiveness in the 
transonic range. This large reduction in amplitude is not evident, how- 
ever, when the stabilizer is used as the manuevering control. 

Each of the computed transfer coefficients shows some irregularity 
at or near the Mach number of 0.92. Otherwise the . following trends are 
evident ; 

The undamped natural frequency con increases with Mach number. 

The danq)lng ratio ^ at 35^000 feet decreases from O.38 at a Ifech num- 
ber of 0.8 to 0.15 at a Mach number of 1.0. The coefficients t, Ki, 
and Ks each show, in general, a decrease with increasing Mach number 
through the transonic speed range. The elevator effectiveness parameter 
^ computed fr^ experimental values of t and Kj. decreases rapidly 
from 0,35 to 0.06 between Mach numbers of 0,8 and 1.0, 

It is shown that the effect of altitude changes on the transfer 
coefficients is such as to cause a decrease in the amplitude and natural 
frequency of the airplane response with increasing altitude. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, California, 
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APPENDIX A 

EELATIONSHIP BETWEEN EXPEEIMENTALLY DETEEMIHED TRAN3EEE 
FUNCTIONS AND THEOEETICAL EqUATIOro OF MOTION 


If the longitudinal equations of motion of the airplane in level 
flight at constant speed are assumed to he of the following form: 

mV(a — q) = ZgO, + Z58 

and 

lyq = Mqq + Maa + I^a + MgB 


they can he solved simultaneously to give 


q _ ^Oq + 

S p® + hp + k 


(Al) 


where p is the operator and where 

du 


k = 


Ma 


(A2) 


mV ly 

^sMa ~ ^6^ 
^q mVI„ 


(A3) 

(A4) 


2.5M5, Ms 

C, = + — 

mVIy ly 


(A5) 


It was shown in this report that the experimentally determined 
pitching— velocity responses may he represented hy equations of the form 


^ - TT 
^ - El 


1 + Tp 


1 + 2C 
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By rearranging terms this equation may be expressed as 

q ^ El + gj_ Tp 
^ p2 + 2C «nP + '^n^ 

which is of the same form as equation (A1) so that 


k = (a6) 

b H 2^o)n (a7) 

C = Kicon^ (a8) 

q. 

= Ki^n^T (a9) 


To obtain the normal— acceleration response the two equations 

of motion may be solved simultaneously for a/B. Then n/B is obtained 
directly from the relation 



and may be expressed as 


n 

B 


p^ + bp + k 


where 


ZoM6 - 


n ml. 




Cs = — i 
n m 
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Inveatigation has shown that for conventional aircraft the term 

is negliglhle and that in most cases has only a s ma ll effect on 

the airplane response over the low frequency range of interest. There- 
fore, n/6 may he expressed simply as 

- = 

6 p^ + hp + k 

This equation is of the same form as the experimentally determined 
transfer function 



and thus is Identical to 

n 
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APPENDIX B 

CALCULATION OF C„ 

“6 

The elevator effectiveness parameter Cjjj^ can be computed directly 

from the experimental results of this report by making use of eq.uations 
developed in appendix A. In equation (A5), the term ZgM^/mVIy is, 
for conventional aircraft, very small as compared to M§/ly, and thus 

M5 may be approximated simply as lyCi • Then from equation (A9) Mr 
can be expressed as Q. ^ 

M5 = ly Ki T 

Finally from the definition of M5, 

_ 2Iy T 

pV^Sc 
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APPENDIX C 

ESTIMATION OF ALTITUDE EFFECTS 


From the definitions of the coefficients obtained from the theoret- 
ical equations of motion, it is possible to estimate the effects of 
changes in altitude on the airplane frequency response. If the Mach 
number remains constant as altitude and density change, then the flight 
velocity will vary directly with the speed of sound (V = Ma). By neg- 
lecting some small terms and by assuming that no changes occur in the 
nondimensional derivatives, these coefficients may be expressed in terms 
of p and a as follows; 


k 


b 



P pV^Sc 
(omitting 


\ 

Sri; 


Zql Mq+I^ 


pvs 

2 




+ C 


“d 



~ Pa 




mVI, 


4ml, 




Cl = 


q I 


y 


21 

5 

(omitting ) 

mVI„ 


Co^ = - ~p=a‘ 


From the identities that define the relations between the various 
coefficients, it can be seen that 


T 1 

pa 
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El ~ Pa 
Zj2 ~ pa^ 
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Figure 1.— Photograph of the test airplane. 
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Figure 2- Two- view drawing of the test airplane. 
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Time, t, sec 

Figure 5— Sample flight records of normal acceleration, pitching velocity, and 
elevator angle at Mach number of 0.8 i . 
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Mach number 


Figure 4.~ Variation of trim elevator angle with Mach number 
and maximum elevator deflection attained during each 
maneuver. 
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Figure 5 -Variation of trim angle of attack with Mach number 
and the maximum and minimum angles following each 
elevator pulse . 
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Figure 6.~ Variation of the first-order term (l+riw) with 
nondimensionat frequency tu . 
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(a) Amplitude . 


Figure 7- The variation of the second -order term 



with nondimensiona! 


frequency — for various values of damping ratio. 
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(b) Phase angle. 


Figure 7 - Concluded . 
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Frequency, u, radians/sec 


(a) Amplitude ratio versus frequency. 

Figure 8. - Airplane pitching-veiocity frequency response to elevator 
control inputs for various flight Mach numbers. 
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Figure 8 -Concluded. 
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Frequency, at, radians/sec 

(a) Amplitude ratio versus frequency. 


Figure 9.- Airplane normal-acceleration frequency response to elevator 
control inputs for various flight Mach numbers. 
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Figure 9- Concluded . 
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Frequency, w, radians/sec 


Figure !0.- Airplane pitching -velocity frequency response to stabilizer 
control inputs for various flight Mach numbers. 
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Figure //.- Variation of the ratio of stabUzer to elevator 

d 

effectiveness, ——r , with Mach number. 

• d,, ’ 
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Figure 12 — Variation of damping ratio and undamped natural 

frequency with Mach number. 
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Figure !3- Variation of the characteristic time,r, and of the 
coefficients K/ and K 2 ^ith Mach number. 
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Figure 14.- Variation of the parameter K^r and elevator effective- 
ness Cm vfith Mach number. 
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Figure !5- Variation of transfer coefficients ,K, ^r,and Kg 
with attitude for constant Mach number and constant lift 
coefficient. 
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Frequency, co, radians /sec 


Figure !6.- Typical pitching-veiacity frequency- respanse 
variation with attitude at Mach number of 0.8! . 
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